Justification of the effective fractional Coulomb energy and the extended 

Brinkman-Rice picture * 
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In order to calculate the effective mass of quasiparticles for a strongly correlated metallic system 
in which the number of carriers, n, is less than that of atoms (or lattices), I, the metallic system 
is averaged by one effective charge per atom over all atomic sites. The effective charge of carriers, 
< e >= (n/l)e — pe, and the on-site effective Coulomb repulsion, U = p 2 (^r) = p 2 kU c , are defined 
and justified by means of measurement, p is band filling and k is the correlation strength. The effec- 
tive mass, = 1 _ a 4 , calculated by the Gutzwiller variational theory, is regarded as the average 
of the true effective mass in the Brinkman-Rice (BR) picture and is the effect of measurement. The 
true effective mass has the same value irrespective of p, and is not measured experimentally. The 
true correlation strength of the BR picture is evaluated as 0.90 < kbr < 1 for Sri-zLa^TiOs and 
0.92 < kbr < 1 for YBCO. High-T c superconductivity for YBCO is attributed to the true effective 
mass caused by the large kbr value. The effective mass explains the metal-insulator transition of 
the first order on band filling. Furthermore, the pseudogap is predicted in this system. 



I. INTRODUCTION 

Ever since a correlation-driven metal-Mott-insulator 
transition 1 of the first order (called the "Mott transi- 
tion") was predicted in 1949, the Mott transition has 
never been completely solved. The Hubbard model 2 , one 
of the theories explaining the Mott transition, showed 
the second-order metal-insulator transition(MIT) by in- 
creasing the on-site Coulomb repulsion U, whereas the 
Brinkman-Rice (BR) picture 3 , based on Gutzwiller's 
variational theory 4 , explained the first-order transition at 
U/U c =l for the first time. The infinite dimcnsional(d = 
oo) Hubbard model suggested that the MIT occurs while 
going from U/W to U c /W, and that the metallic side ex- 
hibits the BR picture. 5 Here W is the band width. The 
d = oo Hubbard model and the BR picture were built 
on a metallic system with an electronic structure of one 
electron per atom. However the metallic system which 
Gutzwiller assumed is a general one without limiting the 
number of electrons per atom . Furthermore, observa- 
tions and current understanding of the MIT driven by 
correlation effects were reviewed. 6,7 Validity of the BR 
picture and the Mott transition was also discussed by 
theorists, which is because the BR picture has not been 
proved experimentally. 8,9 

Experimentally, for the coefficient 7 of the heat ca- 
pacity of Sri-^La^TiOs (SLTO) which is proportional 
to the effective mass, the sharp first-order transition was 
observed between #=0.95 and 1. 10,11 This transition was 
regarded as the Mott transition on band filling, which 
could not be explained by the above mentioned Hub- 
bard models and the BR picture. This is because the 
structure of the metallic system of SLTO is not one elec- 
tron per atom. Hence, when the number of carriers is 
less than the number of atoms, the MIT of first-order on 
band filling has remained a theoretical problem. In par- 
ticular, the metallic property near the Mott transition is 
still controversial. 12-14 Furthermore the effective mass of 



quasiparticles as a function of band filling is necessary to 
take account of high T c superconductivity. 15 

Recently, an insight on both the effective mass and the 
Mott transition was proposed on the basis of U = p 2 (^r) 
and experimental data. 16 However the justification of 
U = /0 2 (^-) and exact calculations for the effective mass 
were not given. Thus the problem of the Mott transition 
on band filling remains unanswered. 

In this paper, when the number of carriers is less than 
the number of atoms, the necessity of the development 
of < e >= pe is given by introducing the instability of 
the charge-density-wave potential energy, and the justi- 
fication of < e >= pe and U = p 2 (^r) is given by means 
of measurement. The effective mass of quasiparticles is 
exactly calculated on the basis of the Gutzwiller varia- 
tional theory. Furthermore, true correlation strengths of 
the BR picture, not given in reference 16, are evaluated 
for Sri-zLa^TiOa and YBCO. 



II. EFFECTIVE FRACTIONAL COULOMB 
ENERGY 

A strongly correlated metallic system of the s band 
structure is assumed for a d — 00 dimensional simple- 
cubic lattice. Let n and / be the number of electrons (or 
carriers) and the number of atoms (or lattices), respec- 
tively. In the case of one electron per atom, i.e., n = I, 
the metallic system is a metal and the existence probabil- 
ity (P = n/l = p= band filling) of electrons on nearest- 
neighbor sites is one. The on-site Coulomb repulsion is 
always given by U = U'={^r) and an electron on an atom 
has a spin, as shown in Fig. 1 (a). However, in the case of 
n < I by doping an element to a base insulator or metal, 
U is determined by probability. The metallic system is 
quite complicated, as shown in Fig. 1 (b). Four types of 
regions for the system can be distinguished as possible ex- 
treme examples. Region A in Fig. 1 (b) has no electrons 



1 



on its atomic sites, which corresponds to a normal insu- 
lator. Region C has the metallic structure of one electron 
per atom. Regions B and D have a charge-density-wave 
(CDW) structure unlike the assumed cubic. This must 
be regarded as a CDW insulator with the CDW-energy 
gap 17 depending on the local CDW-potential energy. 17 ' 18 
Moreover, even when one electron on an atom is removed 
in Fig. 1 (a), both nearest-neighbor sites of the atom 
without an electron and the atom itself also are regarded 
as region B. The CDW-potential energy is defined by 
Vcdw — —E p (Qi — Qj) 2 , where Qi and Qj are charges 
irrespective of spins on i and j sites, respectively. E p 
is the small-polaron binding energy which is a constant 
of electron-phonon coupling. The potential energy is de- 
rived by breathing-mode distortion (or frozen-breathing 
mode) due to a charge disproportionation (6Q = Qi — Qj) 
between nearest-neighbor sites. The CDW-energy gap is 
regarded as a pseudogap (or pinned CDW gap) for this 
system. Here SQ is called the electronic structure factor 
in real-space. SQ=0 indicates that the electronic struc- 
ture is metallic, while SQ ^0 suggests that it is insulat- 
ing. In addition, Vcdw can also be expressed in terms of 
band filling as follows: Vcdw = ~E p (l — p) 2 = —E p p 2 , 
where p = 1 — p = n^/l (or n& = I — n), n& is the number 
of bound charges bounded by Vcdw- 18 As p increases, 
p decreases. This indicates that the pseudogap in the 
CDW insulating phase decreases as the metallic phase 
increases, because the pseudogap depends upon Vcdw- 

Conversely, the pseudogap, as observed by optical 
methods, 19-22 photoemission spectroscopy, 23 and heat 
capacity measurement 24,25 for high-T c superconductors, 
necessarily occurs if the number of carriers is less than the 
number of atoms. Therefore, the metallic system is in- 
homogeneous and cannot be self-consistently represented 
in fc-space, i.e., U and the density of states of the system 
are not given. Thus, in order to overcome this difficulty, 
charges on atoms in the metal phase (region C) have to 
be averaged over sites, as shown in Fig. 1 (c). The on- 
site effective charge is Qi — Qj =< e >= (n/l)e = pe, 
on average over sites. The metallic system, then, is a 
metal, because Vcdw = due to 5Q=0. This system is 
analogous to the case of one electron per atom, such as 
Fig. 1 (a). Therefore, the effective Coulomb energy is 
given by U = p 2 U' , as shown in Fig. 1 (c). 

To illustrate the physical meaning of P = (< e > 
/e) = p < 1, then p « 1, p = 1 — p < 1. For the 
CDW insulating side, p « 1 indicates that there are a 
few doubly occupied atoms such as in region D, while, 
for the metallic side, p < 1 indicates that a number of 
carriers exist. In other words, this suggests that region 
C is extremely wide. 

The insulating and metallic sides correspond to two 
phases, which are attributed to the metal-insulator in- 
stability (or instability of the CDW-potential energy) 18 
at p=l (half filling). This indicates that metals with the 
electronic structure of Fig. 1 (a) are not synthesized. 
Junod et al. 26 suggested that even the best samples are 
not 100% superconducting, which supports the metal- 



insulator instability 18 . Thus synthetic metals composed 
of several atoms always have the electronic structure of 
two phases such as Fig. 1 (b), which is a necessity of the 
development of < e >= pe. 

To justify the fractional charge of carriers of the metal- 
lic system with two phases, the concept of measurement 
is necessary. In the case of measuring charge of carriers 
in the metallic system of Fig. 1 (b) , the measured charge 
becomes the effective fractional one, < e >= pe, with 
the average meaning of the system because experimen- 
tal data for the metallic system are expectation values of 
statistical averages. In the case of not measuring charge, 
the charge of carriers in region C in Fig. 1 (b) remains 
the true elementary one not observed by means of mea- 
surement. 

Accordingly, for the metallic system of the n<l case, 
the effective fractional Coulomb energy, U = p 2 U' , is 
defined by using < e >= pe and justified by means of 
measurement. 

Furthermore, U = p 2 U' , deduced on the basis of the 
CDW concept, can be applied to even cuprate systems 
such as Lai^Sr^C^O^LSCO). This is because the MIT 
with x from the antiferromagnetic insulator, La2Cu04, 
to a metal of LSCO is continuous although there is a 
controversy as for the MIT. 

III. CALCULATION OF THE EFFECTIVE MASS 

For the metallic system regarded as a real synthetic 
metal with an electronic structure such as Fig. 1 (b), the 
effective mass of quasiparticles needs to be calculated. 
Hamiltonians of the metallic system can be considered 
as follows. Hamiltonian, H, is given by 

H = H 1 +H 2 , (1) 
Hi =^( a i T a feT +a{ l aki)e k + U^2 a^a^ag^, (2) 

k g 

H 2 = -J2E P (Q t -Q ] ) 2 , (3) 

where a' k1 . and a' ^ are the creation operators for electrons 
in the Bloch state k and the Wannier state g, respectively, 
and e/c is the kinetic energy when U=0. H\ and H 2 are 
Hamiltonians of the metallic region C and the CDW- 
insulator regions B and D, respectively. In the case of 
Fig. 1 (a) and (c), the Hamiltonian is reduced to H\ 
because H 2 disappears due to 5Q = 0, and the on-site 
Coulomb energy is given by U — p 2 U' . H\ is consistent 
with the Hamiltonian used in the Gutzwiller variational 
theory 4 . 

In order to calculate the effective mass of quasiparti- 
cles and the ground-state energy for a strongly correlated 
metallic system, the Gutzwiller variational theory 4,27-29 
is used. Hi is supposed to describe the metallic system. 
The wave function is written as 
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l*> = ^l*o>. 



(4) 



where I^q) is the wave function when U = 0, v is the 
number of doubly occupied atoms, and < r\ < 1 is 
variation. 4 The expectation value of Hi is regarded to be 



Although the following calculations were performed by 
Brinkman and Rice, 3 the calculations are applied to the 
effective mass. In the case of Fig. f (c), by applying Eq. 
(9) to Eq. (6) and by minimizing it with respect to v, 
the number of the doubly occupied atoms is obtained as 



(H) = 



(5) 



^4 (1+ 8i^4 (1 -L7-)' 



The second part of the equation is simply given by 
Pi}Pii \^)— Uv because |* ) is an eigenstate of 
the number operator pi^Pii- The first part is dealt 
with by assuming that the motion of the up-spin elec- 
trons is essentially independent of the behavior of the 
down-spin particles (and vice versa). By minimization 
with respect to r/, Gutzwiller obtained an extremely sim- 
ple result for the ground-state energy, namely, 

E g /l = <7 T (^,p lT ,p 4i )e T +qi(v,Pn,Pn)e l + Uv. (6) 

Here, 



Sfc<fc F €fc < 



(7) 



is the average energy of the o spins without correlation 
and ek is the kinetic energy in Hi , with the zero of energy 
chosen so that J^k^k—O- ej is equal to ej^. 

The discontinuities, q a , in the single-particle occupa- 
tion number at the Fermi surface are given by 



P)(l -Pa~ P-a + V) + \/{p-a = v)vj 
P<r(l- Pa) 



where p a = |p, < p <1 and p-j- = p^. 27 ' 28 This, calcu- 
lated by Ogawa et. al. 27 who simplified the Gutzwiller 
variational theory, is in the context of the Gutzwiller vari- 
ational theory. Eq. (8) is a function of p a and v irrespec- 
tive of the quantity of charges. This can be analyzed in 
two cases of p=l and < p <1 , because Gutzwiller did 
not limit the number of electrons on an atom for the 
metallic system. 
In the case of p=l, 
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(9) 



This was described in the BR picture. 

In the case of < p < I, two kinds of q a can be 
considered. One is Eq. (8) when the electronic structure 
is SQ ^0, as shown in Fig. I (b). However, Eq. (8) can 
not be applied to the metallic system, as mentioned in the 
above section. The other is Eq. (9) when the electronic 
structure is SQ=0, as shown in Fig. I (c). Eq. (9) is 
obtained from substituting p a in Eq. (8) with p' a . Here, 
p' = (n'/l) = 1 and p' a = i, because the number of the 
effective charges, n' , is equal to I. It should be noted that 
the metallic system with less than one electron per atom, 
as shown in Fig. I (c), is mathematically consistent with 
that with one electron per atom, as shown in Fig. 1 (a). 
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where L7 c =8|e| because of e = ef + <0, U = p 2 U' 
and U 1 = kII c . < k<1 is the correlation strength. By 
applying Eq. (10) to Eq. (9) again, the effective mass is 
given by 



Qa = 
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(11) 



Although the separate conditions are < p<l and 
< k<1, m* is defined under the combined condition 
< up 2 < 1 and is an average of the true effective mass 
in the BR picture for metal phase (region C in Fig. 1 
(b)). The effective mass increases as it approaches k=1 
and p=l. For k^O and p^O, the effective mass decreases 
and, finally, the metallic system undergoes a normal (or 
band-type) metal-insulator transition; this transition is 
continus. The system at np 2 = 1 can be regarded as the 
insulating state which is the paramagnetic insulator be- 
cause v — 0. At a k value (not one), the MIT from a 
(^)metal at a p value of just below p=l to the insulator at 
both p=l and k=1 is the first-order transition on band 
filling. This has been called the Mott transition by a 
lot of scientists including author. However, this is theo- 
retically not the Mott transition which is an first-order 
transition from a value of U to U c in a metal with the elec- 
tronic structure of one electron per atom at p=l, as given 
in the BR picture. The Mott transition does not occur 
in real crystals because a perfect single-phase metal with 
p=l is not made. Conversely, by hole doping (or elec- 
tron doping to a metallic system with hole carriers) of 
a very low concentration, the first-order transition from 
the insulator with D=0 to a metal can be interpreted as 
an abrupt breakdown of the balanced critical Coulomb 
interaction, U c , between electrons. Then, the U c value in 
the insulator reduces to a U value in a metal phase and 
an insulating phase produces due to doping of opposite 
charges, as shown in Fig. 1(d). This first-order tran- 
sition with band filling is very important result found 
in this picture, which differs from the continuous (Mott- 
Hubbard or second-order) transition by a large U given 
by the Hubbard theory. 

In order to obtain the expectation value of the energy 
in the (paramagnetic) ground state, Eqs. (10) and (11) 
are applied to Eq. (6). E g is given by 
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E g /l = e(l - K P 2 ) 



(12) 
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As U/U c — up 2 approaches one, E g goes to zero. 

In addition, the spin susceptibility in the BR picture 
is replaced by 



/4iV(0) 



(1 - iiV(0)K/9 2 t/ ( 



1+V i : 

(l + K^)* > 



(13) 



where N(0) is the density of states at the Fermi sur- 
face and hb is the Bohr magneton. The susceptibility is 
proportional to the effective mass which allows the en- 
hancement of Xs- 

In the above picture, to*, E g /l and \s are the averages 
of true effective values in region C in Fig. 1 (b) which 
is described by the BR picture, and are justified only by 
means of measurement. The true effective ones in the BR 
picture are not measured experimentally, as mentioned 
in an above section. In particular, the magnitude of the 
true effective mass in the BR picture has the same value 
regardless of the extent of region C, while the measured 
effective mass depends upon the extent of region C. 



IV. CONCLUSION 

In order to evaluate true correlation strengths of 
Sri-^LazTiOa and YBCO, explanations not given in the 
previous paper 16 are appended. The experimental data 
of the heat capacity of Sri^La^TiOa was well fitted by 
Eq. (11) with k=1 for all x below x = p=0.95. K=l 
indicates that the effective mass has a constant value 
although p changes. The increase of x to x—0.95 corre- 
sponds to the increase of region C. Because region C is de- 
scribed by the BR picture, the true correlation strength, 
kbr = U/U c , of the BR picture can be found. When it 
is assumed that the extent of the metal phase (region C) 
at x—0.95 is closely the same as that at p—1, a value of 
m*/m = 1/(1 — p 4 ) at x = p=0.95 and k—1 is approxi- 
mately equal to that of m*/m = 1/(1 — n 2 BR ) in the BR 
picture. Then kbr = (0.95) 2 =0.90 is obtained, which 
indicates that Lai-^Sr^TiOs is very strongly correlated. 
Moreover, the decrease of the effective mass from ir=:0.95 
to x—0 is the effect of measurement not the true effect. 

In the case of YBCOg+p at p—0.96, as shown in Fig. 1 
(c) and (d) of reference 16, the true correlation strength 
is determined as kbr=0-92 by the same above method. 
High-T c superconductivity for YBCO is attributed to the 
true effective mass caused by the large kbr value. The 
heat capacity data measured by Loram et al. 24 seems to 
be explained by the above picture, too. 

As another experimental result, the first-order transi- 
tion at £=0.02 for /i-BaNb^Tii-^Oa was observed, which 
may also be in the context of the above picture. 30 The 
validity of the BR picture is indirectly proved through 
the above picture. Furthermore, the combination of the 
above picture and the BR picture 3 is called the extended 



I would like to acknowledge Dr. Kwang-Yong Kang for 
providing the research environment necessary for this re- 
search. 



This was merged in a full-paper presented at NATO Ad- 
vanced Research Workshop on New Trends in Supercon- 
ductivity held at Yalta in Uk raine at Sept. 16-20 of 2001. 
Please see :ond-mat/0110112. 
kimht45@hotmail.com. 
htkim@etri.re.kr.. 

N. F. Mott, Metal-Insulator Transitions, Chapter 3, (Tay- 
lor & Frances, 2nd edition, 1990). 

1 J. Hubbard, (III) ibid 281, (1965) 401. 

1 W. F. Brinkman and T. M. Rice, Phys. Rev. B2, (1970) 
4302. 

1 M. C. Gutzwiller, Phys. Rev. 137, (1965) A1726. 

' X. Y. Zhang, M. J. Rozenberg, and G. Kotliar, Phys. Rev. 

Lett. 70, (1993) 1666. 
' A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 

Rev. Mod. Phys. 68, (1996) 13. 

M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, (1998) 1039. 
! T. Moriya, BUTSURI (edited by Physical Society of 

Japan), Vol 54, No. 1, (1999) 48 (Japanese). 
' H. Fukuyama, M. Imada, and T. Moriya, BUTSURI 

(edited by Physical Society of Japan), Vol 54, No. 2, (1999) 

123( Japanese). 

1 Y. Tokura, Y. Taguchi, Y. Okada, Y. Fujishima, T. Arima, 
K. Kumagai, and Y. Iye, Phys. Rev. Lett. 70, (1993) 2126. 
K. Kumagai, T. Suzuki, Y. Taguchi, Y. Okada, Y. Fu- 
jishima, and Y. Tokura, Phys. Rev. B 48, (1993) 7636. 

! K. Morikawa, T. Mizokawa, K. Kobayashi, A. Fujimori, H. 
Eisaki, S. Uchida, F. Iga, and Y. Nishihara, Phys. Rev. B 
52, (1995) 13711. 

1 1. H. Inoue, O. Goto, H. Makino, N. E. Hussey, and M. 
Ishikawa, Phys. Rev. B 58, (1998) 4372. 

1 H. Makino, I. H. Inoue, M. J. Rozenberg, I. Hase, Y. Aiura, 
and S. Onari, Phys. Rev. B 58, (1998) 4384. 

' C. C. Tsuei, C. C. Chi, D. M. Newns, P. C. Pattnaik, and 
M. Daumling, Phys. Rev. Lett. 69, (1992) 2134. 

; Hyun-Tak K im, Physica C341-348, (2000) 259: |cond- 



mat/000100£ 



BR picture defined in < up < 1. 

ACKNOWLEDGEMENTS 



T. M. Rice and L. Sneddon, Phys. Rev. Lett. 47, (1981) 
689. 

Hyun-Tak Kim, Phys. Rev. B 54, (1996) 90. 
T. Timusk and B. Statt, Rep. Prog. Phys. 62, (1999) 61. 
S. H. Blanton, R. T. Collins, K. H. Kelleher, L. D. Rotter, 
Z. Schlesinger, D. G. Hinks, and Y. Zheng, Phys. Rev. B 
47, (1993) 996. 

M. A. Karlow, S. L. Cooper, A. L. Kotz, M. V. Kelvin, P. 
D. Han, and D. A. Payne, Phys. Rev. B 48, (1993) 6499. 
Hyun-Tak Kim, H. Uwe, and H. Minami, Advances in Su- 
perconductivity VI (Springer- Verlag, Tokyo, 1994), P. 191. 
H. Ding, T. Yokaya, J. C. Campuzano, T. Takahasi, M. 
Randeria, M. R. Norman, T. Mochiku, K. Kadowaki, and 



4 



J. Giapinzakis, Nature 382, (1996) 51. 

24 J. W. Loram, K. A. Mirza, J. R. Cooper, W. Y. Liang, and 
J. M. Wade, Journal of Superconductivity 7, (1994) 243. 

25 J. W. Loram, J. L. Luo, J. R. Cooper, W. Y. Liang, and 
J. L. Tallon, Physica C341-348, (2000) 831. 

26 A. Junod, D. Eckert, T. Graf, G. Triscone, and J. Muller, 
Physica C 162-164, (1989) 482. 

27 Tohr Ogawa, Kunihiko Kanda, and Takeo Matsubara, 
Prog. Theor. Phys. 53, (1975) 614. 

28 Dieter Vollhardt, Rev. Mod. Phys. 56, (1984) 99. 

29 Patrick Fazekas, Lecture Notes on Electron Correlation and 
Magnetism, Chapter 9, (World Scientific Co., 1999). 

30 M. N. Khan, Hyun-Tak Kim, H. Minami, and H. Uwc, 
Materials Letters 47, (2001) 95. 



(a) n=l 
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(c) n<l • charge <e> = (nil) e 




U=«exe>/r> =(n/l) 2 <ee/r>=p 2 U' 



FIG. 1. (a) In the case of one electron per atom, the on-site 
Coulomb repulsion is given by U — { — ) — U' . (b) In the 
n < I case, the four possible electronic structures are region 
A (insulator), region C (metal), and regions B and D (CDW 
insulator). Here, n corresponds to the number of carriers in 
region C. (c) The metallic structure in the case of less than 
one electron per atom is shown. U = (n//) 2 (^-) = p 2 U'. 
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